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Abstract 

The zero-point energy of a conducting spherical shell is evaluated by imposing 
boundary conditions on the potential ^4^, and on the ghost fields. The scheme 

X" 

^ • requires that temporal and tangential components of perturbations should 

vanish at the boundary, jointly with the gauge-averaging functional, first cho- 
sen of the Lorenz type. Gauge invariance of such boundary conditions is then 
obtained provided that the ghost fields vanish at the boundary. Normal and 
longitudinal modes of the potential obey an entangled system of eigenvalue 
equations, whose solution is a linear combination of Bessel functions under 
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the above assumptions, and with the help of the Feynman choice for a dimen- 
sionless gauge parameter. Interestingly, ghost modes cancel exactly the con- 
tribution to the Casimir energy resulting from transverse and temporal modes 
of jointly with the decoupled normal mode of A^. Moreover, normal and 
longitudinal components of Af^ for the interior and the exterior problem give 
a result in complete agreement with the one first found by Boyer, who studied 
instead boundary conditions involving TE and TM modes of the electromag- 
netic field. The coupled eigenvalue equations for perturbative modes of the 
potential are also analyzed in the axial gauge, and for arbitrary values of the 
gauge parameter. The set of modes which contribute to the Casimir energy 
is then drastically changed, and comparison with the case of a flat boundary 
sheds some light on the key features of the Casimir energy in non-covariant 
gauges. 
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Typeset using REVT^ 



1. INTRODUCTION 



One of the most striking applications of quantum field theory to physical processes 
consists of the analysis of zero-point energies. Indeed, it is by now well known that 
zero-point energies of quantized fields are infinite, whereas suitable differences in zero-point 
energies are finite and measurable. A good example is provided by the quantum theory of 
the electromagnetic field. In that case, the energy of two perfectly conducting plates of area 
A and separation d is found to be^ 

720^3 ^ ^ 

One thus has an attraction between such plates. Further to this effect, the zero-point energy 
of a conducting spherical shell was calculated in Ref. 2 for the first time, giving, however, a 
different sign with respect to the case of parallel plates: 

A£;(it:)^o.o9^, (1.2) 

where R is the radius of the sphere. The result (1.2) has been confirmed, independently, 
by all authors working in this research field. A powerful tool turns out to be the method 
relying on the use of the generalized ^-function. The boundary conditions of the problem 
(see below) are of the Dirichlet or Robin type. The corresponding Laplace-type operators, 
say Q, are then self-adjoint, and the analytic continuation of their ^-function: 

CQ(s)=Tri.(Q-^) (1.3) 

may be used to obtain the value of AE{R).^'^ We here omit all details concerning ^-function 
calculations, and we focus instead on the problem of boundary conditions. 

Let r, 9, (p be the spherical coordinates, which are appropriate for the analysis of a con- 
ducting spherical shell. The boundary conditions require that the tangential components of 
the electric field should vanish, i.e., 

[E0]dM = 0, (1.4) 
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[E^]dM = 0. (1.5) 

These boundary conditions are gauge- invariant, since they involve some components of the 
electromagnetic-field tensor: 

F^, = W^A, - W^A^, (1.6) 

which is gauge-invariant whenever the connection V is torsion-free, as we always assume. 
The modes of the electromagnetic field are then split into transverse electric (TE) or mag- 
netic multipole:^ 

r-E = Er = 0, (1.7) 

and transverse magnetic (TM) or electric multipole:^ 

r-B = Br = 0. (1.8) 

On using the notation in Ref. 9 for spherical Bessel functions, and the formulae therein for 
Eg and E^, one finds that, in the TE case, Eqs. (1.4) and (1.5) lead to 

Aai{kr) + Bini{kr)]^^ = Q, (1.9) 

while in the TM case the vanishing of Eg and E^ at the boundary implies that 
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^^yr{Ciji{kr) + Dini{kr))) 



= 0. (1.10) 

dM 



Of course, if the background includes the point r = 0, the coefficients Bi and Di should 
be set to zero for all values of I to obtain a regular solution. Such a singularity is instead 
avoided if one studies the annular region in between two concentric spheres. 

The boundary conditions in the form (1.9) and (1-10) are sufficient to implement the 
algorithm of the ^-function. However, they do not tell us anything about the contribution 
of the ghost field. Indeed, if one follows a path-integral approach to the quantization of 
Maxwell theory, it is well-known that the second-order operator acting on perturbations 
when the Lorenz gauge- averaging functional is chosen turns out to be, in a fiat background,'^ 



P^.u = -g^^uO + (i - ^) V^V,. (1.11) 

With a standard notation, g is the background metric, a is a dimensionless parameter, □ is 
the D'Alembert operator 

where A is the Laplace operator (in our problem, A is considered on a disk). This split of the 
□ operator leads eventually to the eigenvalue equations for the Laplace operator acting on 
the temporal, normal and tangential components of the potential, after using the equations 
(2.1)~(2.3) of Sec. 2. At this stage, the potential, with its gauge transformations 

^A^ = A^ + V^e, (1.12) 

is viewed as a more fundamental object. Some care, however, is then necessary to ensure 
gauge invariance of the whole set of boundary conditions. For example, if one imposes the 
boundary conditions 

[At]aM = 0, (1.13) 
[Ae]9M = 0, (1.14) 

K]9M = 0, (1.15) 

this is enough to ensure that the boundary conditions (1.4) and (1.5) hold. However, within 
this framework, one has to impose yet another condition. By virtue of Eq. (1.12), the 
desired boundary condition involves the gauge function: 

[e]dM = 0. (1.16) 

Equation (1.16) ensures that the boundary conditions (1.13)-(1.15) are preserved under the 
gauge transformations (1.12). What happens is that e is expanded in harmonics on the 
two-sphere according to the relation 
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oo I 

e{t,r,9,ip)=J2 E eiir)Yimi9,ip)e' 

1=0 m=-l 

After a gauge transformation, one finds from (1.12) 



iuit 



[1.17) 



''At\dM — [AtldM = [Vjej^M, 



;i.i8) 



(1.19) 



(1.20) 



and by virtue of (1.17) one has 



[1.21) 



dM 



oo I 

J2 J2 ei{r)Yim,9i^,ip)e 

1=0 m=-l 



J dM 



[1.22) 



dM 



oo I 



J2 J2 ei{r)Yim,^{e,ip)e 



iuit 



=0 m= 



dM 



[1.23) 



Thus, if ei{r) is set to zero at the boundary VZ, the right-hand sides of (1.21)-(1.23) vanish 
at dM. But 

[ei{r)]9M = 0Wl 

is precisely the condition which ensures the vanishing of e at dM, and the proof of our 
statement is completed. 

At this stage, the only boundary condition on Ar whose preservation under the trans- 
formation (1.12) is again guaranteed by Eq. (1.16) is the vanishing of the gauge-averaging 
functional, say ^{A), at the boundary:^ 



MA)]aM = 0. 



:i.24) 



On choosing the Lorenz term ^l{A) = V^Ai,, Eq. (1.24) leads to 
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d_ 
dr 




J dM 



0. 



(1.25) 



It is only upon considering the joint effect of Eqs. (1.13)-(1.16), (1.24) and (1.25) tliat 
tlie wliole set of boundary conditions becomes gauge- invariant. Tliis scheme is also BRST- 
invariant (cf. Ref. 10). More precisely, Eq. (1.16) should be replaced, in the quantum theory, 
by a boundary condition on both the real and the imaginary part of a complex- valued ghost 
field, say'' 



Even this terminology is a bit loose, since ghost fields are really independent of each other, 
as stressed by DeWitt.^^ In the calculations of the following sections, it will be enough to 
consider a real-valued gauge function obeying Eq. (1.16) at the boundary, and then multiply 
the resulting contribution to the zero-point energy by —2, bearing in mind Eq. (1.26) and the 
fermionic nature of ghost fields. A similar analysis of boundary conditions for the Casimir 
effect can be found in Refs. 12 and 13. Note also that our boundary conditions on the 
potential lead to the fulfillment of all "perfect conductor" boundary conditions studied in 
Ref. 2, which involve tangential components of the electric field and the normal component 
of the magnetic field. 

The plan of our paper is as follows. Section 2 derives the form of basis functions for the 
electromagnetic potential. Section 3 presents the eigenvalue conditions for the perturbation 
modes occurring in the Casimir calculation in the presence of spherical symmetry, and 
describes how the Casimir energy of a conducting spherical shell can be obtained from such 
eigenvalue conditions. The eigenvalue equations in the axial gauge are studied in detail in 
Sec. 4. Results and open problems are described in Sec. 5, relevant details are given in the 
Appendices. 



HsM = W\aM = 0. 



(1.26) 
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2. BASIS FUNCTIONS FOR THE ELECTROMAGNETIC POTENTIAL 



We know from Sec. 1 that, in our problem, the electromagnetic potential has a temporal 
component At, a normal component A^, and tangential components, say A^ (hereafter, k 
refers to 9 and ip). They can all be expanded in harmonics on the two-sphere, according to 
the standard formulas 

oo I 

A{t,r,0,<f)^Yl E ai{r)Yim{9,^)e^', (2.1) 

1=0 m=-l 



oo I 

Ar{t,r,9,^) = E E k{r)YU9,^)e'-\ (2.2) 

1=0 m=-l 



oo I 

Ak{t,r,9,ip) = Y, E [Q(r)afcll^(^,<^) + ci/(^)^fep5^>/m(^,<^)]e*"*. (2.3) 

1=1 m=—l 

This means that we are performing a Fourier analysis of the components of the electromag- 
netic potential. The two terms in square brackets of Eq. (2.3) refer to longitudinal and 
transverse modes, respectively. On setting a = 1 in Eq. (1.11), the evaluation of basis 
functions for electromagnetic perturbations can be performed after studying the action of 
the operator —g^u^ on the components (2.1)-(2.3). For this purpose, we perform the ana- 
lytic continuation uj — > iM, which makes it possible to express the basis functions in terms 
of modified Bessel functions (of course, one could work equally well with cu, which leads 
instead to ordinary Bessel functions). This can be shown by using the formulas (A2)-(A4) 
of Appendix A when d — 2. For the temporal component, one deals with the Laplacian 
acting on a scalar field on the three-dimensional disk: 

('"A^), = ^ + H^ + 1(«A^)„ (2.4) 

which leads to the eigenvalue equation 

2d 1(1 + 1 



+ 



dr'^ r dr r 



2 



ai = M^ai. (2.5) 
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Hereafter, we omit for simplicity any subscript for M^, and the notation for the modes will 
make it sufficiently clear which spectrum is studied. The solution of Eq. (2.5) which is 
regular at r = is thus found to be 

1 



ai{r) = -j=Ii+y2{Mr), 



(2.6) 



up to an unessential multiplicative constant. 

The action of ^^^A on the component Ar normal to the two-sphere is (from Eq. (A3)) 



(2.7) 



where the stroke | denotes two-dimensional covariant differentiation on a two-sphere of unit 
radius. Last, the Laplacian on tangential components takes the form (cf. Eq. (A4)) 



gj.2 



+ AA)k-\Ak + -dkAr 



(2.8) 



By virtue of the expansions (2.2) and (2.3), jointly with the properties (A5) and (A6) of 
spherical harmonics, Eqs. (2.7) and (2.8) lead to the eigenvalue equation 



_ /(/ + !) 



di = M^di 



(2.9) 



for transverse modes, jointly with entangled eigenvalue equations for normal and longitudinal 
modes (here / > 1): 



^ ^ 2_d _ (1(1 + 1) + 2) 
dr"^ r dr 



21(1 + 1) , 
hi + ^ , ' ci = M\, 



(2.10) 



_ 1(1 + 1] 



Q + -hi = M'ci. 
r 



(2.11) 



The mode 6o(^) is instead decoupled, and is proportional to '^■v^^'') ^]-^g interior problem. 
It is indeed well known that gauge modes of the Maxwell field obey a coupled set of eigenvalue 
equations. In arbitrary gauges, one cannot decouple these modes. This can be proved by 
trying to put in diagonal form the 2x2 operator matrix acting on the modes hi and q (cf. 
Refs. 7 and 15). In our problem, however, with our choice of gauge-averaging functional and 



gauge parameter, gauge modes can be disentangled, and a simpler method to achieve this 
exists. For this purpose, we point out that, since the background is flat, if gauge modes 
can be decoupled, they can only reduce to linear combinations of Bessel functions, i.e., 

, , , BJMr) 

bi{r) = (2-12) 

and 

Q(r) = C(i/)S,(Mr)v^. (2.13) 

With our notation, C(z/) is some constant depending on u, which is obtained in turn from 
To find u and C(z^), we insert the ansatz (2.12) and (2.13) into the system of equations 
(2.10) and (2.11), and we require that the resulting equations should be of Bessel type for 
B^{Mr), i.e.. 



This leads to two algebraic equations for u 



.2. 



jy^ ^^ + 1(1 + 1) -2l{l + l)C, (2.15) 



and 



+ 1(1 + 1)-^. (2.16) 
By comparison, one thus finds a second-order algebraic equation for C: 

l{l + l)C^ -C -1^0, (2.17) 



whose roots are 



C+ = y, (2.18) 



and 



"--wry 
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The corresponding values of v are 



(2.20) 



and 



(2.21) 



Hence one finds the basis functions for normal and longitudinal perturbations in the interior 
problem in the form 



Ii+3/2{Mr) , Ii-i/2{Mr) 
bi[r) = ai^i '— h a2,i- 



whereas, from Eq. (2.9), transverse modes read 

di{r) = 7,+i/2 (Mr) v^. 



(2.22) 



(2.23) 



(2.24) 



Last, but not least, ghost modes obey an eigenvalue equation analogous to (2.5) (see com- 
ments after (1.26)), i.e. 

2d i{i + iy 



dr^ r dr 



(2.25) 



and hence they read 



1 



ei{r) = —Ii+i/2{Mr). 



(2.26) 



In the exterior problem, i.e., for r greater than the two-sphere radius R, one has simply 
to replace the modified Bessel functions of first kind in Eqs. (2.6), (2.22)-(2.24) and (2.26) 
by modified Bessel functions of second kind, to ensure regularity at infinity. 
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3. EIGENVALUE CONDITIONS AND CASIMIR ENERGY 



In our problem, ghost modes are of course decoupled from the modes for the electromag- 
netic potential which occur in the expansions (2.1)-(2.3). Nevertheless, this does not mean 
that they do not play a role in the Casimir energy calculation. By contrast, we already 
know from Sec. 1 (see comments after (1.16)) that boundary conditions on the ghost are 
strictly necessary to ensure gauge invariance of the boundary conditions on the potential. It 
is then clear that such boundary conditions, combined with the differential equation (2.25), 
lead to a ghost spectrum whose contribution to the Casimir energy can only be obtained 
after a detailed calculation (e.g. Green- function approach, or C- function regularization) . As 
is pointed out in Ref. 12, it is not surprising that ghost terms are important, since one has 
already included effects of other degrees of freedom which should be compensated for. This 
issue is further clarified by the analysis of the original Casimir problem: two perfectly con- 
ducting parallel plates. In a covariant formalism, one has to consider the energy-momentum 
tensor for ghosts, which is found to give a non-vanishing contribution to the average energy 
density in vacuum (see Sec. 4.5 of Ref. 7). After taking into account boundary conditions 
for the potential entirely analogous to our Eqs. (1.13)-(1.15) and (1-24), one then finds 
a renormalized value of the zero-point energy in complete agreement with the result first 
found by Casimir.^ 

By virtue of Eq. (1.13), the modes a/(r) obey homogeneous Dirichlet conditions: 

[ai{r)]eM = Oyi>0. (3.1) 
Moreover, Eq. (1.25) implies that the modes bi obey the boundary conditions 

= 0,V/>0. (3.2) 

dM 

Last, the modes q and rf/, being the tangential modes, obey Dirichlet boundary conditions 
(cf. Eqs. (1.14) and (1.15)) 

[ci{r)]9M = [di{r)]9M = 0, > 1. (3.3) 
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On taking into account how the modes are expressed in terms of Bessel functions (see Eqs. 
(2.6), (2.22)-(2.24)), one thus finds five sets of eigenvalue conditions for the interior and 
exterior problems, respectively: 

(i) Temporal modes: 



/;+l/2(Mi?)=0,VZ>0, 



(3.4a) 



Ki+y2{MR)=0,Wl>0. 



(3.4b) 



(ii) Decoupled normal mode: 



d_ 
dr 



= 0, 



r=R 



(3.5a) 



d_ 
dr 



0. 



r=R 



(3.5b) 



(iii) Coupled longitudinal and normal modes (here u = I + 3/2): 

{v - l/2)Il{MR)h_2{MR) + 3(1/ - l)M^/^_2(Mi?) 



MR 



+ 3/2)Il_2{MR)I,{MR) = 0, 



(3.6a) 



[u - l/2)Kl{MR)K,.2{MR) + 3(z/ - 1)^^4^^K,_2{MR) 



+ 3/2)Kl_^{MR)K^{MR) = 0. 



MR 



(3.6b) 



(iv) Transverse modes: 



7,+i/2(Mi?) = 0,VZ>l, 



(3.7a) 



X;+l/2 (Mi?) =0,V/> 1. 



(3.7b) 



Note that, the boundary being a two-sphere, our "transverse photons" are two-dimensional, 
unlike other physical problems where transverse photons are three-dimensional since the 
boundary is a three-sphere (or another three-surface). 
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(v) Ghost modes (multiplying their C- function by -2): 

7i+i/2(Mi?) = 0,VZ>0, 



(3.8a) 



Ki+i/2{MR) =0,WI> 0. (3.8b) 
The eigenvalue conditions (3.5a) and (3.5b) can be re-expressed in the form 

h/2{MR) = 0, (3.9) 

Ki/2{MR) = 0. (3.10) 

It is thus clear that, by construction, the contribution of Eqs. (3.9) and (3.10) to the Casimir 
energy of a conducting spherical shell cancels exactly the joint effect of Eqs. (3.4a), (3.4b), 
(3.7a), (3.7b), (3.8a) and (3.8b), bearing in mind the fermionic nature of ghost fields. In 
general, each set of boundary conditions involving a set of positive eigenvalues, say {A^}, 
contributes to the Casimir energy in a way which is clarified by the ^-function method, 
because the regularized ground-state energy is defined by the equation (for Re(s) > Sq — 2) 

Eo{s) ^ -\ E(AO^-^ = -^C [s - i) (3.11) 

which is later analytically continued to the value s = in the complex-s plane. Here is 
the usual mass parameter and C, is the ^-function of the positive-definite elliptic operator B 
with discrete spectrum {\k}'- 

CB(s) = Tri2(i3-^) = ^(A,)-^ (3.12) 

k 

In other words, the regularized ground-state energy is equal to — |Cb(~1/2)- Although the 
eigenvalues are known only implicitly, the form of the function occurring in the mode-by- 
mode expression of the boundary conditions leads eventually to £'o(0).^'^ 

The non-trivial part of the analysis is represented by Eqs. (3.6a) and (3.6b). These 
are obtained by imposing the Robin boundary conditions for normal modes, and Dirichlet 
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conditions for longitudinal modes. To find nontrivial solutions of the resulting linear systems 
of two equations in the unknowns ai^i and a2^i, the determinants of the matrices of coefficients 
should vanish. This leads to Eqs. (3.6a) and (3.6b). At this stage, it is more convenient 
to re-express such equations in terms of Bessel functions of order I + 1/2. On using the 
standard recurrence relations among Bessel functions and their first derivatives, one thus 
finds the following equivalent forms of eigenvalue conditions: 



n+i/2 



(MR) 



0,V/ > 1, 



(3.13) 



1+1/2 



(MR) 



2MR 



(MR) 



0,V/ > 1. 



(3.14) 



Thus, the contribution of the coupled normal and longitudinal modes splits into the sum of 
contributions of two scalar fields obeying Dirichlet and Robin boundary conditions, respec- 
tively, with the I — mode omitted. This corresponds exactly to the contributions of TE 
and TM modes (Eqs. (1.9) and (1.10)), and gives the same contribution as the one found 
by Boyer.^ 



4. AXIAL GAUGE 

As is well known, on using the Faddeev- Popov formalism, one performs Gaussian averages 
over gauge functionals, say x^, by adding to the original Lagrangian a gauge- averaging term 
X^PnuX^j where is any constant invertible matrix. ^-"^ This has the effect of turning the 
original operator on field perturbations into a new operator which has the advantage of 
being non-degenerate. It can also become of Laplace type for suitable choices of some gauge 
parameters. Moreover, the result of a one-loop calculation is expected to be x~ ^'^^ P~ 
independent, although no rigorous proof exists on manifolds with boundary.^ In particular, 
for Maxwell theory, P^i, reduces to a 1 x 1 matrix, i.e. a real-valued parameter, and x^ 
reduces to the familiar covariant or noncovariant gauges for Maxwell theory, e.g. Lorenz, 
Coulomb, axial, ... . It is therefore quite important to understand the key features of 
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the Casimir-energy calculations in all such gauges. For this purpose, we consider a gauge- 
averaging functional of the axial type, i.e. 

= n^A^, (4.1) 

where ■nF' is the unit normal vector field — (0, 1, 0, 0). The resulting gauge-field operator 
is found to be, in our flat background, 

= -gi'^u + V'^V'' + -n^n\ (4.2) 

OL 

Note that, unlike the case of Lorenz gauge, the a parameter is dimensionful and has dimen- 
sion [length]^ (see Ref. 17). 

Now we impose again the boundary condition according to which the gauge- averaging 
functional should vanish at dM: 

\^{A)\qm = [n''A^]aM = [A]aM = 0, (4.3) 

which imphes that all bi modes vanish at the boundary (cf. Eq. (3.2)). 

A further consequence of the axial gauge may be derived by acting on the field equations 

P^M^ = (4.4) 

with the operations of covariant differentiation and contraction with the unit normal, i.e. 

V^P'^-A, = 0, (4.5) 

n^P'^M, = 0. (4.6) 

Eq. (4.5) leads, in flat space, to the differential equation 

dAr 2 , , , 

—^ + -Ar^ 0. 4.7 
or r 

This flrst-order equation leads to 6/ modes having the form 

b, = ^. (4.8) 
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Thus, by virtue of the boundary condition (4.3), the modes bi vanish everywhere, and hence 
Ar vanishes identically if the axial gauge-averaging functional is chosen with such boundary 
conditions. 

Moreover, Eq. (4.6) leads to the equation 



lUJ 



9-4, 1 a . It 



which implies, upon making the analytic continuation cu — > iM, 

- K[±l) ^ ^ 0. (4.10) 
dr r^ dr 

At this stage, the transverse modes di obey again the eigenvalue equation (2.9), whereas the 
remaining set of modes obey differential equations which, from Eq. (4.4), are found to be 



d^ ^2 d l{l + 1) 



dr^ r dr 



ai - M^-^^^^ci = 0, (4.11) 



Mai - M^ci = 0, (4.12) 



dr'^ 

ao = 0. (4.13) 

In particular, Eq. (4.13) is obtained from Eq. (4.12) (when / = 0), which is a reduced form 
of the Eq. Pki^A'^ — upon bearing in mind that bi modes vanish everywhere. 

Last, but not least, one should consider the ghost operator, which, in the axial gauge, is 
found to be 

g = -|. (4.14) 

This leads to ghost modes having the form 

Si = So,ie-''r (4.15) 

On the other hand, following the arguments developed in Sec. 1, one can prove that, also in 
the axial gauge, the ghost field should vanish at the boundary, to ensure gauge invariance 
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of the whole set of boundary conditions on the potential. It is then clear, from Eq. (4.15), 
that ghost modes vanish everywhere in the axial gauge. 

On studying the system (4.10)-(4.12) one has first to prove that these three equations 
are compatible. This is indeed the case, because differentiation with respect to r of Eq. 
(4.10) leads to a second-order equation which, upon expressing ^ from Eq. (4.10) and 
from Eq. (4.12), is found to coincide with Eq. (4.11). Thus, we have a system of two 
second-order differential equations for two functions o; and q. However, these functions 
are not independent, in that they are connected by Eq. (4.10). Hence for every value of 
/ one has one degree of freedom instead of two. Finally, we have for every / two degrees 
of freedom, one resulting from Eqs. (4.10)-(4.12), and another, i.e. the transverse mode 
di- Thus, an estimate of the number of degrees of freedom giving contribution to Casimir 
energy coincides with that in other gauges. Moreover, the parameter a does not affect the 
Casimir energy, since a does not occur in any of the eigenvalue equations. 

Unfortunately, we cannot obtain the exact form of the solutions of Eqs. (4.11) and (4.12) 

in terms of special functions (e.g. Bessel or hypergeometric) . This crucial point can be made 

precise by remarking that, if it were possible to disentangle the system (4.11) and (4.12), 

one could find some functions ai, (3i,Vi,Wi, say, such that the 2x2 matrix 

1 v;\ /Ai BA n ai\ 

Wi 1 I \Ci Dj VA 1 / 
has no non- vanishing off-diagonal elements, where the operators Ai, Bi,Ci, Di are the ones 

occurring in Eqs. (4.11) and (4.12) (cf. Refs. 7,14). For example, the first off-diagonal 

element of such matrix is the operator 

Aiai + Bi + Vi(Ciai + Di). 

On setting to zero the coefficients of and ^ one finds that ai = — V; = where ao,/ 
is a constant. But it is then impossible to set to zero the "potential" term of this operator, 
i.e. its purely multiplicative part. 

Nevertheless, there is some evidence that the axial gauge may be consistently used to 
evaluate the Casimir energy. For this purpose we find it helpful to consider a simpler 
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problem, i.e. the Casimir energy in the axial gauge for the case of flat boundary. In this 
case the basis functions are plane waves 

Af, = ^0,^6^^*^-^+'=^^+'=^^-'^*), (4.16) 

where the admissible values of kx,ky, kz are determined by the boundary conditions, which 
are here taken to be analogous to the case of a curved boundary. Let us choose the conducting 
boundaries parallel to the x- and y-axes, while the vector is directed along the z-axis. 
Then the condition of compatibility of field equations in the axial gauge is reduced to 



f-k'' 


00 kx 


00 ky 


Ojkz \ 


00 kx 


kl + kl- 


kxky 


kxkz 


00 ky 


kx ky 


kl + kl-oj-" 


kykz 


\ojkz 


kxkz 


ky k^ 


kl + kl-uj^ + y 



= 0. (4.17) 



Direct calculation shows that this determinant is equal to 

D ^ --kUou' - k'Y . (4.18) 

Hence we have reproduced the correct dispersion relation between energy ou and wave number 
k (to every admissible value of k there correspond two contributions to the Casimir energy 
of the form oo =\ k \). Of course, no non- vanishing ghost modes exist, once that the axial 
type gauge- averaging functional is set to zero at the boundary (cf. (4.3), (4.14) and (4.15)). 
Interestingly, on imposing the boundary conditions, one can set to zero the Az component 
of the electromagnetic potential as was done with Ar in the spherical case. In this case the 
compatibility condition of field equations is reduced to the vanishing of the determinant of 
a 3 X 3 matrix, obtained by omitting the fourth row and the fourth column in Eq. (4.17). 
One can easily see that the determinant of this 3x3 matrix coincides with that of the 4x4 
matrix up to a multiplicative factor ^, which does not affect the dispersion relation. 

To sum up, a complete correspondence can be established between the key features of 
the axial gauge in the cases of flat and curved boundary: the a parameter does not affect 
the Casimir energy, the component of the potential orthogonal to the boundary vanishes, 
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ghost modes vanish, and only two independent degrees of freedom contribute. On going 
from the fiat to the curved case, however, the analysis of the dispersion relation is replaced 
by the problem of finding explicit solutions of Eqs. (4.10)-(4.12), with the corresponding 
eigenvalue conditions. This last technical problem goes beyond the present capabihties of 
the authors, but the exact results and complete correspondences established so far seem to 
add evidence in favour of a complete solution being in sight in the near future. 

5. CONCLUDING REMARKS 

Scientific knowledge makes progress not only when an entirely new result is first obtained, 
but also when alternative derivations of such properties become available, possibly showing 
that further important problems are in sight. The contributions of our paper fall within 
the latter framework. For this purpose, we have studied an approach to the evaluation of 
the zero-point energy of a conducting spherical shell which relies on a careful investigation 
of the potential and of ghost fields, with the corresponding set of boundary conditions 
on Afj, perturbations and ghost modes. When Boyer first developed his calculation,^ the 
formalism of ghost fields for the quantization of gauge fields and gravitation had just been 
developed, and hence it is quite natural that, in the first series of papers on Casimir 
energies, the ghost contribution was not considered, since the emphasis was always put on 
TE and TM modes for the electromagnetic field. On the other hand, the Casimir energy is 
a peculiar property of the quantum theory, and an approach via path-integral quantization 
regards the potential and the ghost as more fundamental.^'^^ This is indeed necessary to 
appreciate that Maxwell theory is a gauge theory. Some of these issues were studied in Refs. 
7,12,13,21,22, including calculations with ghosts in covariant gauges, but, to our knowledge, 
an explicit mode-by-mode analysis of the and ghost contributions in problems with 
spherical symmetry was still lacking in the literature. The contributions of our investigation 
are as follows. 

First, the basis functions for the electromagnetic potential have been found explicitly 
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when the Lorenz gauge-averaging functional is used. The temporal, normal, longitudinal and 
transverse modes have been shown to be linear combinations of spherical Bessel functions. 
Second, it has been proved that transverse modes of the potential are, by themselves, unable 
to reproduce the correct value for the Casimir energy of a conducting spherical shell. Third, 
it is exactly the effect of coupled longitudinal and normal modes of ^4^ which is responsible 
for the value of /\E given in Eq. (1.2), following Boyer^ (see also Refs. 3,4,6). This 
adds evidence in favour of physical degrees of freedom for gauge theories being a concept 
crucially depending on the particular problem under consideration and on the boundary 
conditions. Further study is necessary to understand this property. ^"^'^^ Fourth, ghost modes 
play a nontrivial role as well, in that they cancel the contribution resulting from transverse, 
decoupled and temporal modes of the potential. Fifth, the axial gauge- averaging functional 
has been used to study the Casimir energy for Boyer's problem. Unlike the case of the 
Lorenz gauge, ghost modes and normal modes are found to vanish, and one easily proves 
that the result is independent of the a parameter. A complete comparison with the case 
of flat boundary has also been performed, getting insight into the problems of independent 
degrees of freedom and non- vanishing contributions to the Casimir energy. 

Indeed, recent investigations of Euclidean Maxwell theory in quantum cosmological back- 
grounds had already shown that longitudinal, normal and ghost modes are all essential to 
obtain the value of the conformal anomaly and of the one-loop effective action. ^'^^'^^'^^ In 
other words, the mode-by-mode evaluation of the Faddeev-Popov path integral in the one- 
loop approximation agrees with the result obtained, for the same path integral, by using 
the method of conformal variations with the corresponding Schwinger-DeWitt asymptotic 
expansion. ^'^^ For example, in a portion of flat Euclidean four-space bounded by two con- 
centric three-spheres, the contribution to C(0) of the decoupled normal mode cancels exactly 
the one resulting from transverse modes. Moreover, on the four-dimensional disk, trans- 
verse modes give contributions of opposite signs to C(0) with magnetic or electric boundary 
conditions,^ whereas the (^(0) value is the same for both boundary conditions and equal to 
—31/90.^^'^^ Thus, in problems where the background is four-dimensional and the bound- 
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ary is three-dimensional, the covariant ^-function calculations, including ghost modes, give 
results in disagreement with a calculation relying on boundary conditions for the fields 
and reduction to physical degrees of freedom7'^^'^^'^^ Further evidence of the nontrivial role 
played by ghost modes in curved backgrounds had been obtained, much earlier, in Ref. 25. 
Hence, we find it nontrivial that the ghost formalism gives results in complete agreement 
with Boyer's investigation.^ Note also that, in the case of perfectly conducting parallel plates, 
the ghost contribution cancels the one due to tangential components of the potential (see 
Sec. 4.5 of Ref. 7). This differs from the cancellations found in our paper in the presence of 
spherical symmetry. 

The main open problem is now the explicit proof that the Casimir energy is independent 
of the choice of gauge-averaging functional. This task is as difficult as crucial to obtain a 
thorough understanding of the quantized electromagnetic field in the presence of bound- 
ing surfaces. In general, as shown in Sec. 4, one has then to study entangled eigenvalue 
equations for temporal, normal and longitudinal modes. The general solution is not (ob- 
viously) expressed in terms of well known special functions (cf., for example, Ref. 14). A 
satisfactory understanding of this problem is still lacking, while its solution would be of 
great relevance both for the foundations and for the applications of quantum field theory. 
In covariant gauges, coupled eigenvalue equations with arbitrary gauge parameter also lead 
to severe technical problems, which are described in Appendix B. A further set of nontrivial 
applications lies in the analysis of more involved geometries, where spherical symmetry no 
longer holds, or in the investigation of media which are not perfect conductors. Thus, 
there exists increasing evidence that the study of Casimir energies will continue to play an 
important role in quantum field theory in the years to come, and that a formalism relying 
on potentials and ghost fields is, indeed, crucial on the way towards a better understanding 
of quantized fields. 
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APPENDIX A: 

One can consider, in general, a {d+ 1) -dimensional disk with metric 

g^dr0dr + r'^Qa, r e [0, R] , (Al) 

where fl^ is the metric on the d-dimensional sphere of unit radius. The action of the Laplace 
operator on a scalar field, say 4>, is found to be 

A^-(^ + ^|4«a).. (A2) 

whereas the action of the same operator on normal and tangential components of a one-form 
is 

(A^)^ = + + ^ -|) ^- - ^-^p". (A3) 



With our notation, ^'^^ A =|p'^ is the Laplacian on a d-sphere of unit radius. When d — 2, its 
action on spherical harmonics is such that 

A Yim{e, if) = + l)Yim{e, ip). (A5) 

Another useful formula in the course of deriving Eq. (2.11) from Eq. (2.8) is the one 
according to which^^ 

A dkYim{9, 99) = [-/(/ + 1) + l]dkYim{0, (A6) 
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APPENDIX B: 



A naturally occurring question is whether one can prove explicitly that the Casimir 
energy is independent of the value of the a parameter in the operator (1.11). Indeed, from 
Eqs. (2.4), (2.7), (2.8), jointly with the form of the Lorenz gauge-averaging functional: 



+ 



dt dr r 



one finds that 



-Pt.A''= (— + A) A + ^u; 

\ a \a 



dA,,. 1 4 |ik " A \ 



2 

r 



(Bl) 



(B2) 



- Pr^A"" = (u;^ + A)Ar + iu; 1 - 



1\ dA, 



a) dr 



1--\^{^-^ + -A\' + -a)\, (B3) 
a) dr \ dr f"^ ^ r ^ j 



.V / 2 . / . ^\dAt / l\ d (dA. 



As in Sec. 2, we now perform the analytic continuation uj — > iM. For arbitrary values 
of a, this remains a well defined procedure, but it no longer leads to a system admitting 
Bessel type solutions. By contrast, by virtue of the equations (2.1)-(2.3) and (B.2)-(B.4), 
one finds a coupled system of three ordinary differential equations with variable coefficients. 
On defining the column vector 

( ai\ 



u = 



hi 
\ci/ 



jointly with the nine operators 



dr^ r dr r^ \ a J ' 



(B5) 



(B6) 



(B7) 



(B8) 
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.4.,.Mli-l)-i, (B9) 



A - 1 ^ 2 1 ('^^in^i\\^ (BIO) 
a dr^ ar dr \a J r^^ 



A^^(^-l)'(l±ll^,l'l±ll, (BID 

V aj dr a r"* 

^31 = M (^i - l) , (B12) 

^3..(i-l)-f + li, (B13) 

Vq; / otr ar 

^33 = — - (B14) 
the resulting system reads (with summation over the repeated index p = 1, 2, 3) 

AjpUp = M'^Uj, (B15) 

for all j = 1, 2, 3. Remarkably, the temporal, normal and longitudinal modes are now all 
coupled, but the system of equations is so complicated that no explicit solution in terms of 
well known special functions can be obtained, as far as the authors can see. Nevertheless, if 
q; 7^ 1, regular basis functions still exist, since r = is a Fuchsian singularity of the system 
(B.15). This implies that the desired solutions can be expressed, in the interior problem, in 
the form 

oo 

^ Vi,kr\ 

k=0 

where p is a real-valued parameter, and rji^k are some coefficients. One should also say that, 
in the quantum cosmological backgrounds studied over the last few years,^ the problem is 
less severe, since only normal and longitudinal modes remain coupled for arbitrary values 
of a (but in that case there is no time variable, the background is four-dimensional and the 
boundary is three-dimensional). Only the transverse modes are decoupled, and obey the 
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equation (2.9). Moreover, the modes oq and 60 are coupled as well (unlike the case when 
a — 1), and obey the following system of three eigenvalue equations: 



(P 2d 
dr^ r dr 



a 



ao + M 1 



a 



dr r , 
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a dr"^ ar dr ar"^ 



.a J dr 



.a / \a J dr a r 



(B16) 

(BIT) 
(B18) 



Similar technical problems occur on considering other choices of the gauge-averaging func- 
tional, e.g. 

and hence we only write down the eigenvalue equations (B.15) which are obtained in the 
Lorenz case. 

Since in Sec. 4 we have found it quite helpful to perform a comparison with the analysis of 
a flat boundary, here we describe the same procedure for the Lorenz gauge. On using again 
the plane-wave expansion (4.16) for the electromagnetic potential, one gets the following 

compatibility condition for field equations: 

( up' — — ^uP /3u;kx /3ujky 

(5ujkx k"^ — uj"^ — (5k1 —f3kxky 

(5ujky —f3kxky k"^ —uj"^ — (5ky 

f3uj}z^ l3kx^z i^^y^z 

where 3=1 — -. The direct calculation shows that 

1 



det 



V 



Pcukz \ 
—j3kxkz 
(3kykz 

e-cu^-pkij 



= 0, (B19) 



F = --(u' - ky. 

a 



(B20) 



Hence one can see that the dispersion relation does not depend on a. The same holds for 
ghost modes, where the similar condition reads 



2\2 



{u' - k') 



0. 



(B21) 



Thus, we have obtained again the correct dispersion relation, with the correct power of 
iuj^-k^). 
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